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We study the ground state properties of a superconductor - ferromagnet - insulator trilayer on the 
basis of a Hubbard Model featuring exchange splitting in the ferromagnet and electron - electron 
attraction in the superconductor. We solve the spin - polarized Hartree - Fock - Gorkov equations 
together with the Maxwell's equation (Ampere's law) fully self-consistently. For certain values of the 
exchange splitting we find that a spontaneous spin polarized current is generated in the ground state 
and is intimately related to Andreev bound states at the Fermi level. Moreover, the polarization of 
the current strongly depends on the band filling. 

PACS numbers: 72.25.-b, 74.50.+r, 75.75.+a 
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I. INTRODUCTION 

Recently, the proximity effect between a superconduc- 
tor {SC) and a ferromagnet {FM) has attracted much 
attention because, due to advances in materials growth 
and fabrication techniques well controlled structures, 
in which it can occur became available. Such FM /SC 
hybrid structures are important from the point of view of 
their intrinsic scientific interest, as they allow the study 
of the interplay between ferromagnetism and supercon- 
ductivity 1^ as well as of device applications in such areas 
of technology as magnetoelectronics or quantum com- 
puting g. 

In the present context the proximity, on one hand, 
means the leakage of superconductivity into the ferro- 
magnet, and on the other, the spin polarization of the 
superconductor. In the normal metal - superconductor 
systems the proximity effect has been studied over 30 
years, and is now well understood Q to be governed by 
the Andreev reflection ||^ processes, in which an elec- 
tron at the interface is reflected as a hole while a Cooper 
pair is transfered into superconductor. This mechanism 
makes the movement of electrons and holes in the normal 
metal highly coherent. 

If the normal metal is replaced by a ferromagnet a 
number of new phenomena may occur. For instance one 
may expect a version of the FFLO state, predicted in 
sixties for a bulk superconductor in a strong spin ex- 
change field, by Fulde and Ferrell jl^ and Larkin and 
Ovchinikov [Q, to be realized. In the bulk, as it is well 
known, the exchange field tends to polarize the conduc- 
tion electrons and holes. In particular if these electrons 
are Cooper paired, naively one would expect that either 
the exchange field is too week to break the pairs, or it 
leads to the first order phase transition to the normal 
state. However, as was demonstrated in Ref. [|ll|, for 
certain values of the exchange field a new superconduct- 
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ing depairing state is realized through first order phase 
transition from the BCS state and it transforms continu- 
ously, by a second order phase transition, into the normal 
state as the strength of the exchange field is increased. 
This FFLO state has a spatially dependent order param- 
eter corresponding to the nonzero center of mass motion 
of the Cooper pairs. Another novel feature of this state 
is a current flowing in the ground state. The unpaired 
electrons tend to congregate at one portion of the Fermi 
surface so a quasiparticle current is produced. In order 
to satisfy the Bloch theorem: no current in the ground 
state, a supercurrent, generated by the nonzero value of 
the pairing momentum, flows in opposite direction, and 
the total current is zero. 

Interestingly, similar oscillations of the pairing ampli- 
tude (Cooper pair density) have been predicted 
in a ferromagnet proximity coupled to a superconductor. 
It turns out that these oscillations are responsible for the 
oscillatory dependence of the critical temperature on the 
thickness of the ferromagnetic slab ^ . This effect as 
well as the corresponding oscillations of the density of 
states at the Fermi level ^ , have been observed ex- 
perimentally [p^ [l^ . In this paper we investigate further 
the ramifications of these interesting phenomena. 

One of these is the formation of so called tt junction 
state which is also a fingerprint of the oscillatory be- 
haviour of the pairing amplitude. The 7r-state effect has 
been extensively studied in connection with the high- Tc 
superconductors, where due to spatially inhomogcneity 
of the order parameter, the Cooper pair wave function 
can change sign at the interface and this leads to the for- 
mation of a zero - energy mid - gap state. Remarkably, 
this zero - energy state is unstable to the occurrence of 
a spontaneous current fiowing parallel to and within a 
coherence length of the interface (see eg. and refer- 
ences therein). In what follows we report on our finding 
similar spontaneous currents but of different origin in a 
SC/FM/I trilayer. 

We also find zero energy bound states due to the finite 
size of the FM slab. In fact even for a N/SC hcterostruc- 
ture, in agreement with , we have found states within 
the SC gap. Although the energy of these states depend 
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on the thickness of the normal metal, they never reach 
zero, namely ep. Interestingly, when the normal metal 
is replaced by a ferromagnet, these bound states split, 
and for certain values of the exchange field, they cross 
the Fermi level As we shall show, this circumstance 
leads to a current flow. Clearly, the spontaneous cur- 
rent is strictly related to the zero energy states, as in 
the case of the high - superconductors, but the ori- 
gin of these states is completely different. As mentioned 
above, there are two competing parts to the total current 
in the FFLO state. One is related to the unpaired elec- 
tron movement and the other is a supercurrent. So one 
may expect, that in the case of the FM/SC heterostruc- 
ture, these parts could be spatially unbalanced creating 
regions of net current flow. Indeed this is what we have 
discovered. 

A brief report of our main results has already been 
published |Q. Here we wish to present a more detailed 
and systematic study of the zero energy Andreev bound 
states and the corresponding spontaneous current in the 
FM/SC heterostructures. The paper is organized as fol- 
lows: In Sec. II the simple model which allows for self- 
consistent description of the FM/SC heterostructure is 
introduced. We also derive equations which have to be 
solved self-consistently. In particular, these include the 
Maxwell's equation (Ampere's law) for the coupling of 
the current to the magnetic field. Then we describe some 
technical details concerning the principal layer technique, 
which allows_us to describe a semi-infinite superconduc- 
In Sec. 



tor. 



Ill the nature of the Andreev bound states 



in the ferromagnet is discussed. The spontaneously gen- 
erated current and corresponding magnetic field in the 

In Sec. 



ground state are studied in the Sec. [V 



we 



provide some suggestions of how one might observe the 
spontaneous current and its polarization experimentally. 
Finally, we conclude in Sec. 
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II. 



THE MODEL 



A. Negative U Hubbard Model with exchange 
splitting 

To study spontaneous currents and their polarization 
in the ferromagnet - superconductor heterostructure we 
have adopted a single orbital, nearest neighbour hopping 
Hubbard model with negative U on the SC side and zero 
U otherwise. Additionally we have assumed the site en- 
ergies EiCT to be exchange split in the ferromagnet and 
spin independent in the superconductor. Moreover, the 
simplest geometry allowing for a current flow is the 2D 
system, shown in the Fig.y where the magnetic field in 
one direction, the vector potential and the current in an- 
other and the spatial modulation in a third orthogonal 
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FIG. 1: Schematic view of the (finite thickness) ferromagnet 
- semi-infinite superconductor heterostructure. Directions of 
the magnetic field (B) as well as vector potential (A) and 
current (J) are indicated. 



direction are explicitly indicated. 

In short, our model Hamiltonian is given in the form: 

H = Y^ik] + (Sta - y)5^j]ct^Cj„ + ^ yn„,n,;_<, (1) 



where, in the presence of a vector potential A{r), the 



hopping integral is given by ti 



-te 



-ie /-^ A{r)-dr 



for 



nearest neighbour lattice sites, whose positions are fi 
and r}, and zero otherwise. The site energies £ia are 
on the superconducting side and equal to ^-E'e^cr on the 
ferromagnetic side, /i is the chemical potential, the on- 
site interaction Ui is Us < in the superconductor and 
zero elsewhere, c^, (cio-) are the usual electron creation 
(annihilation) operators and fiic, — c^^Cia is the electron 
number operator. Note that the above description of the 
electrons with charges e includes a coupling to a magnetic 
field B (f) ^ W X A{r) , which is necessary for calculating 
the effects of current on the electronic states. 

In what follows we shall work within Spin - Polar- 
ized - Hartree - Fock - Gorkov (SPHFC) approximation, 
which means that we have approximated the interaction 
term in Hamiltonian (^ by the usual mean field theory 
mapping: n^T^ia ^ {cziCit)cl^cj'^ + {c'l'^cii)ciiCi^. For the 
magnetic field we have chosen the Landau gauge where 
B = {0,0,B^{x)) and hence A = {0,Ay{x),0). Fur- 
thermore, we have assumed that the eflFective SPHFG 
Hamiltonian is periodic in the direction parallel to the 
interface and therefore we work in k space in the y di- 
rection but in real space in the x-direction (see FigJ^). 
Labelling the 'planes' (lines) in Fig.0 by integers n and 
m at each ky point of the Brillouin zone we shall solve 
the following SPHFC layer index (n) matrix equation 
for the retarded Green's function matrix Cm,'m{'-^,ky): 



E 

m'ku 



(2) 



where fg is the Pauli unit matrix and Hnm{ky) is of the 
form: 
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2 ^ex^nm 



A„(5„ 








(3) 



with T± = {tcos{ky ± eAy{n)) + /i)(5„,„ + t(5„,n+i. 

B. Semi - infinite superconductor 

Our system consists of infinite number of layers (n), be- 
cause of fact that we have a semi-infinite superconductor. 



The FM region spreads from n = —Np to 0, while super- 
conductor is defined for n > f (see Fig. In order to de- 
scribe semi-infinite system, we need infinite range matrix 
with block elements M„m(w, ky) = LOToSnm^ Hnmiky) (n, 
m going from —Np through to oo). In short, strictly 
speaking, we have to invert the infinite range matrix: 
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To render the problem traceable, it is useful to define the 
surface Green's function (SGF) |^ as 

G^^U^,ky) = {[M„>^,(^,fc,)]-i}„,„ (5) 

Mn-^f^giuj^ky) being semi-infinite submatrix for n,m > 
Ns (lower right corner of the matrix in Eq.(^). The 
physical meaning of the SGF is that it represents a 
square diagonal subblock corresponding to the semi- 
infinite bulk shown with a shaded region in the Fig. 0. 

Now we can assume, that electronic properties in the 
SG differ from that in the bulk only over the finite dis- 
tance from the interface (1 < n < Ns). Ns can cor- 
respond to the distance of a few SC coherence lengths. 
For n > Ns the electronic properties are that of the bulk 
superconductor. In other words, we assume that the sys- 
tem for n > Ns is homogeneous. Namely, for n > Ns we 
have: 

^ Mn,n{(^, ky) = Mns,Ns ky) 
Mn,n+l{i^,ky) = MMs,Ns + l{^^ky) (6) 

M„ + 1,„(W, fcy) = MMs + l,Ns{^^ky) 



As a consequence the surface Green's function (||) is also 
n-independent and can be obtained from the bulk GF 




G'f{uJ,ky) = GNs.Nsi^i'ky) X 

(1 + Mns+i,Ns{^, ky)GNs.Ns+i{^, ky)y^ C^) 



where G'„,m(w, fcj,) is the Fourier transform of 
the Green's function of the homogeneous super- 
conductor G(uj,kx,ky), namely Gn,m{^,ky) — 
l\t Efe. G(u;,fc:,,/cj,)e-*'=-(^"-^'"). Clearly, the Green's 

function G{Lo,kx,ky) for the case where the vector 
potential is equal to zero is given by the usual result: 
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tj + ^k -A 

-A Lj-Ck 

A2) [ cj + ^k A 

A w-Ck 



(8) 



where ^k = —2t{cos{kx) + cos{ky)) — /i is the single elec- 
tron energy spectrum. 

The homogeneity of the system for n > Ns leads to 



the conclusion, that the inversion of a semi-infinite block 
tridiagonal matrix (^) can be reduced to the inversion of 
the finite tridiagonal matrix: 



M(W, ky) = 
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In the above equation the effect of the semi-infinite ho- 
mogeneous system is contained in the 'embedding poten- 
tial' Fjvsj which is related to the surface Green's function 
through: 

fArs(w,fcy) = MMs,Ns + l{'^iky) X 

G'f{uJ,ky)MNs+l,Ns{^^.ky) (10) 

Evidently, by increasing Ns until Gnmi^^ky) in the 
surface region does not change would constitute an exact 
albeit numerical solution of the Eq. (||). All calculations 
reported here were performed for Ng = 20. However, 
we have checked that larger Ns has no influence on the 
results. 

C. Finite temperature method 

The calculation of the various physical quantities, like 
particle concentration n„, spin polarization m„, order 
parameter A„, etc., needs the evaluation of integrals of 
the product of the Green's function and Fermi distribu- 
tion function as in / dujG{uj)f{uj). While working on the 
real energy axis, we need a large number of integration 
points to get good accuracy, unless there are no singu- 
larities in the density of states. However, in our calcula- 
tions, there is a Van Hove singularity in the 2D DOS as 
well as SCS'-like singularities. In fact we need approx- 
imately 10"^ - 10^ points to get good accuracy, which is 
very time-consuming in self-consistent calculations. The 
situation is similar for summation over the Matsubara 
energies uj^ = {2v + l)TTi/P because of the poles of the 



Fermi function. There is an infinite number of poles, and 
infinite sum does not always converge well, so numerical 
treatment needs a large number of data points. 

To overcome these difficulties we follow Ref. and 
approximate the Fermi distribution function by: 

/'(^) = —TT^N (11) 

(^) +1 

To get the correct behaviour near the Fermi energy 
{Ep = 0) we choose 2N = (3a. One can check, that 
the discrepancy between the approximate function and 
the true one, starts to play role for energies lower than 
2a. If we choose 2a > W, where W = 2t is the half of 
the bandwidth, the approximation works very well. For 
example, for the temperature T = 10~H it is sufficient 
to consider only 150 points to get very good accuracy. 

The poles of the approximate function lie on a circle 
given by 

(^) (") 

The solution of the above equation gives the poles in the 
form: 

In the limit a ^ oo, the poles (|l^) move to the Mat- 
subara energies, and the number of poles, 2N, becomes 
infinite. 
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The calculation of any physical quantity becomes: 



(O) = ± y ReGK)e(2-+i)-/2Ar (^4) 

TT ^ — ' 

where uji, is given by Eg. (p^) . 

D. Self-consistency with the Ampere's law 

As usual in problems of interfaces, a large number 
of equations has to be solved self-consistently. In our 
case the parameters, to be determined self-consistently 
at each layer n, are: the total electron density n„, the 
spin polarization m„, the SC order parameter A„, the 
current parallel to the interface Jy (n) and the vector po- 
tential Ay{n). Most of these quantities can be expressed 
in terms of the corresponding layer - diagonal elements 
of the Green's functions G„„(a;, fc^) (see Eq.(||)). The 
exception is the vector potential, for which we need an 
additional (Maxwell's) equation. These relations are (see 
Eq.(ll)): 

2N-1 



rin = "-«T + "-"i = o X! X! 



ky 1^=0 



Re 



{(Gi^c.,, ky) + Gllico,, fc,))e(2'^+i)-/2^} (15) 



2N-1 



^ ky v=0 



Re 



{(G^U-., ky) - Glliuj., fc,))e(^^+i)-/2^} (16) 

A„ = Un^{Cnl{ky)Cn^{ky)) = 



the end of each self-consistency cycle. We calculated this 
new vector potential by solving numerically Ampere's 

law, ^ji"^^ — —4,'Kjy{x), which for the lattice problem 
at hand, is 

Ay{n + 1) - 2Ay{n) + Ay(n - 1) = -47rJy(n) (19) 

In all our calculations the Eqs.(p^)-(p^) have been 
solved in each iteration step until self-consistency has 
been achieved. The numerical calculations have been per- 
formed at temperature T — lO^^t for 150 energy points, 
120 - ky points and 40 layers (20 ferromagnetic and 20 su- 
perconducting). Typically self-consistency, at the level of 
0.01 % on all densities in Eqs.([l5|)-(|l9|), has been achieved 
after 100 - 200 iterations. 



III. ANDREEV BOUND STATES IN 
FERROMAGNET 



Superconducting and ferromagnetic order 
parameters 



Since we have determined the ferromagnetic (16) and 
superconducting ( |l7|) order parameters on both sides of 
the interface fully self-consistently, we were able to study 
both SC and FM proximity effects. This means, that we 
were able to describe situations in which ferromagnetism 
and superconductivity coexist near the interface. 

Firstly, we want to discuss the magnetic proximity ef- 
fect, namely the entering of the spin polarization into the 
superconductor. The typical behaviour of the spin polar- 
ization (m„ — rt„| — nn\) is plotted in the Fig. |^ as a 
function of the layer index n for a number of exchange 
splittings Ef^x ■ We can see that the spin polarization ex- 



2N-\ 



'^Y.Y. Re{G-(c.„fc,)e(-+i)-/^^} (17) 



ky I'—O 



where Gl^„{uj,ky) is the ij -th element of the matrix 
Green's function, which is the solution of the SPHFG 
equation (Eq. (||)). 

The current, for spin up electrons, in the y-direction 
can be calculated from the relation: 



J, 



aT(i)W = -^^sin{ky -eAy{n)) x 

ky 

Y: Rc{G,\i(3^)(c..,fc,)e(2''+i)W2A'} (18) 



2N-1 



which follows from the continuity equation for the charge 

(e^ = -eK,H]_ = -E, J.,)- 

This current will give rise to a vector potential A„^^ (r) 
which will have to be used to update A{'r) in Eq.(0) at 
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FIG. 2: Magnetization as a function of the distance from the 
interface for a number of exchange fields (Eex) and Us ~ —2, 
which gives As — 0.376 in units of the hopping integral t. 

ponentially decays over the distance of the SC coherence 
lengths ^5, which in the present case is ^5 « 3 in units of 
the lattice constant a. We have checked this by explicit 
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calculations of ^5 for various values of the superconduct- 
ing energy gap. So we can conclude that effect of prox- 
imity of the ferromagnet on the superconductor is very 
similar to case of the Meissner effect, where an external 
magnetic field is excluded from the sample. Note however 
that while the magnetic field is excluded on the spatial 
scale of the penetration depth the effective exchange field 
drops to zero within the distance of the coherence length 
in the superconductor. 

More interesting is the behaviour of the spin polariza- 
tion on the ferromagnetic side of the heterostructure. As 
it is seen from Fig. ^, it oscillates around its bulk value 
with period depending on the exchange splitting E^.^ but 
not on the superconducting gap. However we were not 
able to fit the numerical points to any simple relation 
between parameters of the model. The best fit we have 
is nin (X sin{{1.74:Eex + Q-5)n). Whilst this means that 
the period of the oscillations is a linear function of the 
exchange splitting the above relation is valid only for the 
thickness of the FM slab d — 20 layers. For different 
thicknesses we have different relations and different num- 
bers of period within the slab. For example, for d — 10, 
rUn shows one period less for the corresponding exchange 
splittings Egx- 

We now turn to the superconducting proximity effect, 
namely the leakage of the SC properties into the non- 
superconducting region. In our case, although A„ — 
on the ferromagnetic side, due to the fact that for n < 0, 
Un = 0, the pairing amplitude Xn = {cniCn]) can, in 
general, be non-zero and it usually is. 

The proximity effect is well understood in the 
case of non-magnetic metal - superconductor interface 
(NM/SC) ^ in terms oftheAndreev reflections ||]. The 
Andreev process describes a situation in which an elec- 
tron impinging onto the NM /SC interface is reflected 
back as a hole (with opposite spin) and a Cooper pair is 
created in the SC. Moreover, the movement of the elec- 
trons and holes in the NM close to the interface (over 
the distance of the SC coherence length ^5) is highly cor- 
related. The pairing amplitude, in this case, enters the 
normal metal and decays exponentially over the ^s- 

When the normal metal is replaced by a ferromagnet 
one might naively expect that the proximity effect will be 
suppressed due to the pair breaking effect of the exchange 
held . This suggests that the pairing amplitude Xn 
will decay over a distance much shorter than the super- 
conducting coherence length. Remarkably, as we can see 
in the Fig. H, this is not the case. Instead, we can observe 
a very long range proximity effect with an oscillating pair- 
ing amplitude. This effect has been first noted by A. I. 
Buzdin et al. B and have been studied in several re- 
cent publications j]^, ^ |l5[ . According to them, 
it can be attributed to a FFLO - like phenomena 
in FM/SC heterostructures. The period of the oscilla- 
tions of the Xn depends on the exchange splitting E^x 
(see Fig. ^ and, similarly to the oscillations of m„, does 
not depend on the superconducting energy gap. Clearly 
this behaviour is different from NM / SC case, where the 
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FIG. 3: The pairing amplitude — Xn vs layer index n for three 
values of the exchange field Eex- Note that we plot — Xn which 
corresponds to positive A„ for f/„ > 0. Inset: Comparison 
of the numerical results with the analytical formula Xn oc 
sin[n / £^f) / {n / ^f) for the exchange field E^x = 1.88, which 
gives = 1.06. In fact we had to change by 5% in 
order to get better agreement. The solid line is our fit: X" = 
0.021sm(n/|F - 2.5)/(n/|F), where |f = Cf/1.05. 



proximity effect depends on the SC coherence length ^5. 
Here, it depends only on the properties of the ferromag- 
net. Moreover, the numerically obtained behaviour of 
the pairing amplitude Xn is consistent with the analyti- 
cal formula 

Xn oc sin{n/£,F)/{n/^F) (20) 

with S^p — 2t/Eex being ferromagnetic coherence length. 
Reassuringly this result is fully consistent with those 
of Ref. |l5[ The comparison our numeri- 

cal results with the formula (|2^) for the exchange split- 
ting Eex = 1.88, which gives £,f — 1-06 is depicted 
in the inset of the Fig. ||. In fact we had to change 
£,F by 5%, in order to get better fit. So the solid line 
in the inset of the Fig. corresponds to the formula 
Xn = 0.021sm(n/|F - 2.5)/(n/^F) with |f = Cf/1-05. 
It turns out that as exchange splitting becomes smaller 
the fitting to the analytical formula ( pO| ) becomes better. 
The small discrepancy in the coherence length may come 
from the fact, that we have used bulk value of the ex- 
change field while calculating the coherence length. How- 
ever in our system, the spin polarization oscillates around 
its bulk value, and if we calculate the average value of 
the spin polarization rh — J^^_^_l X]n<o ^ 
tie bit smaller value than the bulk one. So, we conclude, 
that the effective exchange splitting E^x also has to be 
smaller than its bulk value. Indeed, when we used E^x 
to calculate the coherence length, the agreement is im- 
proved. Another reason may be the approximate form of 
the formula for the coherence length in the lattice case 

^i? W 2t/Eex- 



7 



B. Density of states and the Andreev levels 

The proximity effect manifests itself not only in the 
pairing amplitude and the spin polarization characteris- 
tics but also in the density of states {DOS). The layer 
resolved DOS is defined as 

-- J2 ^MGliiUJ + lO+, ky) + G^iu, + Z0+, ky)) (21) 



with Gll^iyU) -I- iO"*", ky) being the ij-th element of the ma- 
trix Green's function - solution of the Eq. (|^). We shall 
also use the density of states integrated over the total 
sample: ptot{^) = HnPn{^)- 

At this point we want to stress that, when we calculate 
thermodynamic quantities (Eq.(|l5|)-(|l8|)) we work on the 
complex energy plane, but calculating the DOS we solve 
Eq. I on the real energy axis only once. In fact we have 
added a small imaginary part (0"'" = 0.005) for numerical 
purposes. The quantities: n„, to„, A„, Ay(n) as deter- 
mined by the previous self-consistent procedure allow us 
to find the Green function on the real energy axis in only 
one step. We also check if the new quantities like n„, m„, 
etc. agree with the old ones. For all such calculations the 
agreement has been found to be excellent. 

While the general solution of Eq. (||) allows for a current 
flow, in this section we show the results for a solution for 
which the current has been constrained to be zero. 

To appreciate the effect of spin polarization in the FM 
slab, let us recall again the case of no exchange splitting. 
If the non-magnetic metal (NM) is in the contact with 
superconductor (SC), the density of states on the NM 
side shows features of the SC DOS, i.e. the pseudo-gap 
opens up around the Fermi energy in the DOS. However 
it never turns into a true gap, namely the NM DOS al- 
ways possess low energy excitations. Moreover, the gap 
in the SC DOS also starts to fill up, due to the proximity 
of the normal metal. So in this case we deal with gap- 
less superconductivity close to interface of the NM / SC 
heterostructure. 

Compared to the above summary, the ferromagnet - 
superconductor interface is enriched with a number of 
new features. For instance, the SC DOS splits, if the 
exchange field is small. For strong enough exchange field, 
SC gap can be suppressed, indicating transition to the 
non-SC (normal) phase. An example of the layer re- 
solved density of states for spin up (solid line) and down 
(dashed) is shown in the Fig. |^. Evidently there is a 
pseudogap in the FM DOS, indicating that supercon- 
ducting correlations are present in the FM layers. This 
effect is particularly strong for the layer index n — 
and —1. On the other hand, the proximity of the fer- 
romagnet leads to existence of the low energy states in 
the SC DOS. This effect is also observed in NM/SC 
heterostructure but here we also notice differences in the 
spin up and spin down DOS on the SC layers close to 




FIG. 4; The layer resolved density of states for the spin 
up (solid line) and spin down (dashed line) electrons for the 
exchange splitting E^x = 0.44. The layer index is shown in 
the panels. 



the interface. These are clear manifestations of the FM 
proximity effect. 

The most interesting physics of the heterostructures 
we consider, is the formation of Andreev bound states 
in the NM or FM layers. They are 'particle in a box' 
like states which can be associated with the semiclassical 
orbits bouncing back and forth between the SC and / 
regions, as depicted in the Fig. ^. As is well known 
in the case of NM layers the reflections at the SC are 
Andreev ones, while those at the / are normal reflections 
|]l9| . The arrows in the Fig. ^ point along the momenta 
of the particle segment of the orbit. Its y component is 
ky which labels the state. Clearly, the collection of states 
labelled in this way forms an Andreev band. 

For the normal metal - superconductor heterostruc- 
ture, Andreev bound states (bands) are symmetrically 
placed with respect to the Fermi energy. Of course the 
position of these states changes with the thickness of the 
sample. As we increase the size of the sample, the An- 
dreev bound states approach the Fermi energy, but they 
will never reach it. In other words, there is no possibility 
of crossing the Fermi level. Somewhat surprisingly, the 
situation is quite different for the FM /SC structure [ p^ . 
In this case the bound state energies can cross the Fermi 
level. In short they can become the zero - energy mid - 



gap states 
in the Fig. 



This situation is illustrated, for k^ 



0, 



where the bound state energies are plotted 
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I FM SC 



FIG. 5: Schematic view of the semi-classical paths of imping- 
ing electrons and reflected holes corresponding to the different 
Andreev bound states. 



against the exchange sphtting. 




FIG. 6: The position of the Andreev bound states vs ex- 
change splitting Eex for the 20 layers thick ferromagnet. The 
solid (dashed) line corresponds to the spin up (down) elec- 
trons. Note that energy is now measured in units of the bulk 
superconductor order parameter As. The basic physics rep- 
resented by this plot is the same as the lu/ As vs d (thickness 
of the FM sample) curves of Ref. |l5| . 

If wc take into account states for different ky (Andreev 
band) , the situation is more complicated as in general the 
states with different symmetry of the wave function can 
mix, and so we are not able to separate Andreev bands so 
easily. The example of such Andreev bands rather than 
states is shown in the Fig. 0, where the density of states 
of the spin up electrons is plotted for a few values of the 
exchange splitting. 

As we shall report presently the fact that the Andreev 
bound states can cross the Fermi level has remarkable 
consequences: it leads to the generation of the sponta- 
neous current in the ground state. 

In fact the situation is rather similar to that of ex- 
otic superconductors (with non-s-wave order parameter) , 
where the surface DOS shows the zero-energy state, lead- 
ing to the spontaneously generated current. However the 
origin of this state is quite different, as there it comes 



0.35 




0.05 I ' ' ' ' 

-0.2 -0.1 0.1 0.2 

CO 



FIG. 7: The total density of states (ptoti i^^) ~ X^„PnT('^) 
of the spin up electrons for various values of the exchange 
splitting Eex indicated on the picture. Note the shift of the 
Andreev bound states (bands) with exchange splitting. The 
DOS for spin down electrons can be easily obtained by reflec- 
tion of the spin up electrons DOS with respect to the Fermi 
energy ep — 0. 



from the symmetry related sign change of the order pa- 
rameter at the surface |Q. Here, in the FM/SC het- 
erostructure, the zero-energy state comes from the os- 
cillatory behaviour of the pairing amplitude (x) in FM , 
which changes its phase by tt each time crosses zero. Nev- 
ertheless it leads to the same effect: a spontaneously gen- 
erated current. 



IV. CURRENT IN THE GROUND STATE 



A. Splitting of the zero energy state 

The most remarkable feature of our calculations is that 
the solution of the SPHFG equations frequently con- 
verges to a solution with the finite current jy{n) even 
though the external vector potential is zero. The gen- 
eration of the spontaneous current is strictly related to 
the crossing of the Andreev bound states through the 
Fermi energy. In the Fig. |^, we have plotted the DOS 
for different values of the exchange splitting for calcula- 
tions where the current was allowed to develop in the y 
- direction. For the E^x — 0.13 we observe that spon- 
taneous current is generated in the ground state. We 
have checked that the solution with the current flowing 
is the true ground state, as it has energy lower than so- 
lution without the current. Physically, in the presence 
of the current, the zero energy states split because of 
the momentum dependence oi a, p ■ v contribution to the 
quasiparticle energies where the velocity vector is defined 
by the current J = env. The example of such a Doppler 
splitting in the DOS is shown in Fig. ^ for the exchange 
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field Eex = 0.13. We can see that, besides the zero energy 
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0.0 I ' ' ' ' ' 1 

-0.4 -0.2 0.2 0.4 

FIG. 8: The Doppler splitting of the zero energy state in the 
density of states ptot^jj) = I]„(pnT('^) + Pni(t^)) caused by 
the spontaneous current for E^x = 0.13. The solid (dashed) 
line corresponds to the solution without (with) the current in 
the ground state. 

state, the other Andreev states also are split. It turns out 
that the splitting of these states is only weakly sensitive 
to the exchange field E^x , but it strongly depends on the 
thickness of the ferromagnetic slab and it decreases as 
the thickness increases. 

Interestingly, we have found a simple relation between 
the splitting of the zero energy state 5 and the vector 
potential Ay{n). This is given by: 

5 « 2etAy (22) 

where Ay is the value of the vector potential Ay{n) av- 
eraged over the FM side of the system only, i.e. Ay = 

Corresponding to the splitting of the zero energy 
states, when the current flows, we find differences in the 
DOS at the Fermi level. In fact, the difference is quite 
small. The main behaviour remains the same. The com- 
parison of the density of states at the Fermi energy for 
both cases (J = and J 7^ 0) is shown in the Fig. ||. 
The DOS for the NM/SC structure (dashed fine with 
circles) is also drawn. As we can see, the main features 
remain more or less the same, in particular the period 
of the oscillations. Similarly like the spin polarization 
(Fig. ^) and pairing amplitude (Fig. ||), the density of 
states at the Fermi level shows decreasing period of the 
oscillations as the exchange splitting increases. 



B. Spontaneous current and the magnetic flux 

As we already discussed, the zero energy state is re- 
sponsible for the generation of the spontaneous current 
in the system. The typical example of such a current, 
flowing parallel to the FM/SC interface, {ji°*{n) = 
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FIG. 9: The layer resolved density of states at the Fermi 
energy ef = for different exchange splittings. The filled 
(empty) symbols correspond to the solution with (without) 
the current. 
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FIG. 10: The total (spontaneous) current jl°*{'n) — jy] {n) + 
jyl{n) flowing parallel to the FM/SC interface for a number 
of exchange splittings. 

jylin) + jyii^)) is shown in the Fig. |l^ for a few val- 
ues of the exchange splitting. Evidently, the behaviour 
of the current, as a function of the layer index, is very 
similar to the density of states at the Fermi level. For ex- 
ample, if we compare the dashed curve with the triangles 
in Fig. |l^ to the dashed or solid line with triangles in 
Fig. 0, we readily notice the similarity in the oscillating 
nature of the current and DOS respectively. 

Another important issue is the distribution of the cur- 
rent through the whole trilayer structure. We find that it 
flows mostly in the positive y direction on ferromagnetic 
side and in the negative direction in the superconductor. 
Notably the total current, integrated over the whole sam- 
ple, is equal to zero within numerical accuracy. This is 
as it should be for the true ground state and found to be 
in the FFLO state, where the current associated with 
the unpaired electrons is balanced by the supercurrent 
flowing in the opposite direction. 

While the total current in the system is equal to 
zero, the current integrated over FM side only {Jpl^ — 
J2n<o Jy'^i^)) a flnite value, almost independent of 
the exchange splitting. But that current is very sensi- 
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tive to the thickness of the FM slab, and for example 
for Np = 10 jy^lj « 0.08 (0.008 per layer), while for 
Nf = 20 we have J^^lj « 0.035 (0.00175 per layer). This 
suggests that effect associated with the spontaneous cur- 
rent is very important in the samples of the small thick- 
nesses and it seems to play minor role if the size of FM 
part of the system is large. 

Obviously, the above current distribution should lead 
to the generation of the magnetic flux through the sam- 
ple. In fact we have found that again, as the current, 
the magnetic flux bounded by the layers only weakly de- 
pends on the exchange splitting but it does change with 
the distance from the interface. For Np — 10 the mag- 
netic flux $*°* = YlnKo is nearly equal to the half of 
the flux quantum <i>o, « 0.45$o- This gives a flux 
per plaquette, associated with the layers n and n -I- 1, 
$„ « 0.045$o- On the other hand, for Np = 20 we have 
« 0.25$o and « 0.0125*o respectively. 

The typical layer dependence of the spontaneous mag- 
netic flux $„ (magnetic flux through a plaquette) is 
shown in the Fig. O. We can see that while there are 




FIG. 11: The (spontaneous) magnetic flux per plaquette as- 
sociated with layer n and n + l for a number of exchange split- 
tings. Each curve corresponds to the related current shown 
in the Fig. |lo[ 

small differences in the behaviour on the FM and SC 
sides of the sample, related to the different current dis- 
tributions (see the Fig. ^0|), the total magnetic flux 
is roughly the same and more or less independent of the 
exchange splitting. 



C. Quasiparticle current vs supercurrent. 

As we have recalled above, in the bulk FFLO state 
the total current vanishes. Moreover, the quasiparti- 
cle current J'^ is exactly equal to the supercurrent J*^ 
(related to the G'^^(w,fcy) in Eq. (||)) with opposite 
sign. This means that in a bulk layered superconduc- 
tor, these two currents cancel each other layer by layer, 
i.e. J'^P{n) = — J'"^(n). Evidently this is the case because 



the exchange field and superconducting pairing potential 
are present in the whole system. By contrast, in the case 
of the FM/SC hcterostructure, the exchange field and 
the pairing potential are spatially separated. In this case 
J'^P{n) is not equal to — J'"^(n) within each layer, but it 
is when integrated over the whole sample, i. e. 



^J^^'(n) = -5]j-(n) 



(23) 



Thus the lack of exact cancellation of the current layer by 
layer leads to a finite current on both sides of the inter- 
face, but zero in the whole sample. An example of such 
behaviour is shown in Fig. 11^. We see, that the quasi- 




FIG. 12: Layered resolved quasiparticle current (solid line) 
and supercurrent (dashed line) for exchange splitting Eex = 
1.88. This corresponds to the total current (solid line with 
circles) shown in Fig. |l^. 

particle part of the total current (solid line) is almost 
in anti-phase to the supercurrent (dashed line). More- 
over, on the superconducting side, the current is carried 
mainly by quasiparticles. This is rather surprising since 
usually one expects current to be related to the Cooper 
pairs. However due to the proximity of the ferromagnet, 
there are also 'normal' particles in the superconductor 
near the interface. These certainly contribute to the to- 
tal current, like in the bulk FFLO state. In our system, 
only the Cooper pairs close to the interface can feel the 
exchange field of the ferromagnet and have moving cen- 
ter of mass. Note that magnetization goes to zero very 
quickly as we go away from the interface (see Fig. ||). 
Clearly this makes the supercurrent equal to zero except 
within a few SC layers. 

Interestingly we also found that whenever the pairing 
amplitude Xn changes sign, the local quasiparticle cur- 
rent as well as the local supercurrent becomes large. In 
other words, the maxima of the J''P{n) and J'''^{n) cor- 
respond to the zeros of the pairing amplitude x„. 



D. The pairing amplitude at I/FM interface 

As we have mentioned above, the ground state does 
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not carry current for all values of the exchange split- 
ting. We have shown, that whenever there is a current 
flowing, there are the Andreev bound states crossing the 
Fermi level. Also we have found another quantity which 
is related to the current flowing. Namely, the pairing am- 
plitude at the FM/I interface (i.e. n = —Np - see Fig. 
|l|). If X-Np changes its sign, the spontaneous current is 
generated. This happens exactly for the same values of 
the exchange splitting, for which a band of the Andreev 
bound states crosses the Fermi energy. Moreover, the 
presence or absence of the current has a dramatic effect 
on the pairing amplitude X-Np ■ This can be seen in the 
Fig. |l^, where X-Np is plotted as a function of the di- 
mensionless parameter <d = 2.79dEex/'n't. As implied by 




FIG. 13: Pairing amplitude at the FM surface opposite to 
superconductor (n — —Np) as a function of the dimensionless 
parameter O — 2.7'ddEex /i^t. 

Eq.|l^, Xn is the real part of a complex order parameter 
which has an amplitude Xn and a phase (pn- Clearly, Fig. 

can be red as Cp-Np flipping from ip'^pf^ where X-Np 
has the same sign as x in 5C, to tp-Np — tt, where X-Np 
has the opposite sign. 

Evidently, the pairing amplitude is pinned to zero by 
the spontaneous current in the system. Introducing the 
dimensionless quantity Q we can say, that system is 
in a TT-like phase (for odd Q) or in the 0-phase (even 
0), in close analogy to studies on the 7r-effect in the 
SC/FM/SC structure by Chtchelkatchev et al. In 
Ref. we have shown a similar picture for the thick- 
ness of the FM slab equal to 10 layers. Here we show 
the corresponding results for Np — 20 layers. To be pre- 
cise in the present case has slightly different value with 
prefactor 2.79 (not 3 as it was for 10 layers case). This 
discrepancy may come from the fact, that effective ex- 
change splitting is little bit sm aller t han its bulk value 
(see discussion in the paragraph III A ). 



ductor, we can expect the spontaneous current to be 
polarized. However, up to now, we have presented re- 
sults for the special case of the particle - hole symmetry 
{ue — 1 or fj, — 0) and we found that the spontaneous 
current was, in fact, not spin - polarized. The reason is 
as follows. Within linear response theory, the total cur- 
rent can be divided into two parts: a diamagnetic one - 
giving the response of the bulk density, and a paramag- 
netic one, which is proportional to the density of states 
at the Fermi energy, as it comes from the deformation of 
the wave function at the Fermi surface |^. Now if we 
note, that for the particle - hole symmetric case (rig — 1), 
the spin up ptotii^p) and spin down ptotiisp) DOS at 
the Fermi energy are equal, we are not surprised that the 
polarization of the current is equal to zero. As we go 
away from the half filling, the difference between spin up 
and spin down DOS {Aptotisp) = Ptot\{£F) - Ptotiisp) 
starts to play a role as it leads to the polarization of the 
current. A typical example of the polarization of the cur- 
rent AJy{n) = Jyiin) — Jyi{n) is shown in the Fig. |l^ 
for concentration of electrons n„ = 0.781 and 0.622. To 
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E. Band filling and the polarization of the current 

As our system consists of ferromagnet and supercon- 



FIG. 14: Spin polarized current AJy{n) — Jy^(n) — Jyi(n) 
as a function of the layer index n for exchange splitting Bex = 
0.44 and concentration of the electrons rie = 0.781 (solid line) 
and 0.622 (dashed line). 

make the two calculations comparable we have adjusted 
the parameters so the value of the SC order parameter 
would be the same. This was achieved by shifting of the 
coupling constant Us to —2.129 for Ue = 0.781, —2.345 
for Tie — 0.622 as well as chemical potential to fi — —0.5 
and —1.0 respectively. 

The total polarization of the current AJ*°* = 
Jj,(n) is of order lO^'^ and does not depend much 
on the concentration of electrons rig. Moreover, we have 
found a very interesting property of this current, namely, 
the sign of the polarization of the current depends on 
whether rie < 1 or > 1. This means, that for = 1 + x 
we get the same results as ior He — 1 — x, except of the 
sign of the polarization of the current. This effect can be 
explained very easily if we recall that polarization of the 
current is proportional to the difference in the DOS at 
the Fermi energy, and one can check, that Aptot{£p) is 
an antisymmetric function under the shift 1 — x ^ 1 + x, 
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i.e. Aptotiep, I - x) = -Aptotisp, 1 + x). 



V. HOW TO OBSERVE THE CURRENT AND 
ITS POLARIZATION ? 

Clearly, the important issue raised by the above results 
is "how to observe experimentally the predicted current 
?" But first of all let's estimate how big this current is. If 
we assume a bandwidth equal to 1 eV, the average cur- 
rent per layer is of order of tenths of pA. In general the 
polarization of the ferromagnet can be large, however it 
cannot be equal to 1. There must be a small number of 
the minority electrons. In other words it doesn't matter 
how big the ratio ^f/Cs we should always observe the 
current. The next point is the thickness of the FM re- 
gion. Namely, the current is bigger for a thinner sample, 
but even for thickness of order of several hundreds of ^i?, 
we expect the current to be easily detectable. 

Clearly, the presence of the current should manifest it- 
self in a number of experiments. For example, since there 
is a net magnetic flux associated with the current, it could 
be detected by SQUID experiments. Another possibil- 
ity is to measure the conductance across the FM j SC 
junction psf . In such experiments the splitting of the 
zero energy state in the DOS is expected to be observed 
and, as indicated by Eq.^ the energy distance between 
the split peaks measures the current flowing. Finally one 
can also use some local probes like STM techniques to 
measure the layer resolved DOS. 

Confirming that the current is polarized is more trou- 
blesome. As we have shown, to get polarized current the 
system has to be away from the rig = 1 point. In our 
case, to avoid a Fermi energy mismatch, both the FM 
and SC are moved from rig = 1. Fortunately, this is only 
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necessary for the ferromagnet. The superconductor can 
be in a particle - hole symmetric state. In short, realistic 
estimations of the spin polarization of the spontaneous 
current require more material specific calculations. 



VI. CONCLUSIONS 

In conclusion, we have studied properties of the insula- 
tor (or vacuum) - ferromagnet - superconductor trilayer. 
We have shown, that such a structure supports Andreev 
bound states forming Andreev bands, whose position can 
be tuned by thickness of the sample or exchange splitting. 
Moreover, when a band crosses the Fermi energy, a spon- 
taneous current (and magnetic flux) is generated in the 
ground state. We have found a relation between the pair- 
ing amplitude in the ferromagnet, which has oscillatory 
behaviour, and the Andreev bound states. Namely, when 
the Andreev band crosses the Fermi energy the pairing 
amplitude at the surface opposite to the FM /SC in- 
terface changes sign, and as long as current flows it is 
pinned to zero. The polarization of the current strongly 
depends on the band flUing and is related to the differ- 
ence in the spin up and spin down density of states at 
the Fermi level. We also have discussed the question of 
the supercurrent vs quasiparticle current and gave some 
experimental clues on how to observe that current. The 
presented results can be, and we have argued that they 
are, attributed to the FFLO effect influenced Andreev 
bound states and can be verifled experimentally. 
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